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^ . Abstract 

JlLi' Significant progress in understanding quarkonium production requires improved 

^ ■ polarization measurements, fully considering the intrinsic multidimensionality 

of the problem. We propose a frame-invariant formalism which minimizes the 
dependence of the measured result on the experimental acceptance, facilitates the 
>• ' comparison with theoretical calculations, and provides a much needed control 

J^ ■ over systematic effects due to detector limitations and analysis biases. This 

QQ ! formalism is a direct and generic consequence of the rotational invariance of the 

CN I dilepton decay distribution and is independent of any assumptions specific to 

lO ' particular models of quarkonium production. 
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Quarkonium polarization measurements should provide key information for the 
understanding of quantum chromodynamics (QCD) pQ, with the competing mecha- 
nisms dominating in the different theoretical approaches to quarkonium production 
leading to very different predictions. In particular, non-relativistic QCD calcula- 
tions [2], dominated by colour-octet components, predict that, at Tevatron energies 
and asymptotically high p^, the directly produced J/tp and ip' mesons should be 
almost fully transversely polarized (angular momentum component Jz = ±1) with 
respect to their own momentum direction (helicity frame), while NLO calculations 
of colour- singlet quarkonium production [3] indicate a strong longitudinal (J^ = 0) 
polarization component. 

The present experimental knowledge on quarkonium polarization is also contra- 
dictory and puzzling. The slightly longitudinal prompt-J/?/' polarization measured 
by CDF |1] along the helicity axis {J/ip momentum direction in the center-of-mass 
of the collision system) is in clear disagreement with the expectations of all existing 
models. The fact that these recent measurements disagree with the results previously 
published by the same experiment |5] adds further confusion to the picture. Equally 
disturbing is the lack of continuity between fixed-target and collider results [6] . Bot- 
tomonium polarization should be easier to interpret theoretically, given that the heav- 
ier bottom quark mass should satisfy the non-relativistic approximation much better 
than the cc case. However, the T(IS') Tevatron measurements [7], available in the 
helicity frame and extending to Pt values around 15 GeV/c, are contradictory: CDF 
indicates unpolarized production; DO indicates longitudinal polarization. The dis- 
crepancy between the two results cannot be justified by their different rapidity win- 
dows. At lower energy and pt, E866 j8] has shown yet a different polarization pattern: 
the T{2S) and T{3S) states have maximal transverse polarization with respect to 
the direction of motion of the colliding hadrons (Collins-Soper frame [9]), with no 
significant dependence on transverse or longitudinal momentum, while the T(IS') is 
only weakly polarized, indicating a dominant role of feed-down contributions. 

To clarify this intricate situation, improved measurements are needed. So far, 
most experiments have presented results based on a fraction of the physical informa- 
tion derivable from the data: only one polarization frame is used and only the polar 
projection of the decay angular distribution is studied. These incomplete results 
prevent model-independent physical conclusions. Moreover, such partial descriptions 
of the observed physical processes reduce the chances of detecting possible biases 
induced by insufficiently controlled systematic effects. The forthcoming LHC mea- 
surements, in particular, would benefit from an improved formalism. In Ref. [6] we 
emphasized the importance of approaching the polarization measurement as a mul- 
tidimensional problem, determining the full angular distribution in more than one 
frame. In this letter we show the existence of a frame-independent relation among 
the observable parameters of the dilepton decay angular distribution. The deter- 
mination of invariant quantities facilitates the comparison between measurements, 
and with theory, reducing the kinematic dependence of the results to their intrinsic 
(and physically relevant) component. Invariant relations can also be used to perform 
self-consistency checks which can expose unaccounted systematic effects or eventual 
biases in the experimental analyses, a crucial advantage given the challenging nature 



of quarkonium polarization measurements due, in particular, to the difficult subtrac- 
tion of the spurious kinematic correlations induced by the detector acceptance. 

We begin with preliminary considerations on the kinematics of the dilepton decay 
of inclusively produced vector mesons. The most general observable distribution is 

W{cos'&,ip) oc -— — — (l + A^cos^^ 

+ A^ sin^ -d cos 2lp + A^^ sin 2'd cos ip) , 

where ■& and ip are the (polar and azimuthal) angles formed by the positive lepton 
with, respectively, the polarization axis z and the production plane xz (containing the 
directions of the colliding particles and of the decaying meson). All experimentally 
definable polarization axes belong to the production plane. A transformation from 
one observation frame (A) to another (B) is a rotation around the y axis hj a 6 angle, 
the parameters changing as: 

(B) ^ Xj) - 3A (B) ^ X^ + A 

'' 1 + A ' ^ 1 + A ' 



X(B) 



Affcos2^-|(A(^^-Ai.^))sin2^ (2) 

^^ 1 + A 

with A = 1 (Ai^) - A(^)) sin^ 6 - \ A^^^ sin 25 . 

Since the magnitude of the "polar anisotropy", A^, never exceeds 1 in any frame, we 
deduce the frame-independent inequalities 

|A^|<^(1 + A^), |A^J<^(1-A^), (3) 

which imply the bounds |A^| < 1 and |A^^| < 1. More interestingly, we can see that 
|A^| < 0.5 when A^ = and must vanish when A^ — )■ —1. 

In general, the transformation of the polarization parameters explicitly depends 
on the production kinematics. Considering, for example, the CoUins-Soper (CS) and 
helicity (HX) frames, the angular terms in Eq. |2] are 

• 2 c- • 2 c- Pt: ^ 

sm ohx->cs — sm ocs->hx — -^tt — 5— — tt , 

sm 2dux-.cs = - sm 2()cs^hx = r,, ^ , — jT ' 

p'^[m'^ +Pt) 

where m, E, p, px and pl are, respectively, the mass, the energy and the total, trans- 
verse and longitudinal momenta of the meson in the center-of-mass of the collision. As 
a result, the observed quarkonium polarization has, in general, an "extrinsic", frame- 
related, kinematic dependence, superimposed on the "intrinsic" physical dependence 
due to the properties of the production processes and their varying mixture. 

Such an extrinsic dependence can introduce artificial differences between the po- 
larization results obtained by experiments probing different acceptance windows. Fig- 
ure [1] shows how natural T polarizations A^ = +1 and —1 in the CS frame (with 
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Figure 1: Kinematic dependence of the T(IS') decay angular distribution seen in 
the HX frame, for natural polarizations A^ = +l (a-b-c) and A^ = — 1 (d-e-f) in the 
CS frame. The curves correspond to different rapidity intervals; from the solid line: 
\y\ < 0.6 (CDF), \y\ < 0.9 (ALICE), \y\ < 1.8 (DO), \y\ < 2.5 (ATLAS and CMS), 
2 < |y| < 5 (LHCb). For simplicity the event populations were generated flat in 
rapidity. The sign of X^^p depends on the definition of the y axis of the polarization 

frame, here taken as sign(pL)(-f"i x P'2)/\P'i x P'2\, where -P'1,2 are the momenta of 
the colliding particles in the meson's rest frame. 



A(^ = A^(^ = and no intrinsic kinematic dependence) translate into different p^- 
dependent polarizations measured in the HX frame in different rapidity acceptances. 
In this simple case, a common choice of the "natural" frame (CS) by all experiments 
would avoid such a misleading differentiation of results. In general, however, it may 
be impossible to find one frame providing a simple representation of the quarkonium 
polarization scenario. This is shown in Fig. |2l where we consider, for illustration, 
that 60% of the T events have natural polarization A^ = +1 in the CS frame and 
the remaining fraction has A^ = +1 in the HX frame. While the polarizations of the 
two event subsamples are intrinsically independent of the production kinematics, in 
neither frame will measurements performed in different transverse and longitudinal 
momenta windows find identical results. 
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Figure 2: Kinematic dependence of the T(IS') decay angular distribution seen in the 
HX (a-b-c) and CS (d-e-f ) frames, when 60% (40%) of the events have full transverse 
polarization in the CS (HX) frame. The curves represent measurements in different 
acceptance ranges, as detailed in Fig. [TJ 



Also the comparison between experimental data and theory must consider the 
experimental acceptance and efficiency. Experiments measure the net polarization 
of the specific cocktail of quarkonium events accepted by the detector, trigger and 
analysis cuts. If the polarization depends on the kinematics, the measured angu- 
lar parameters depend on the effective population of collected events in the probed 
phase space window. Two experiments covering the same kinematic interval may 
find different average polarizations if they have different acceptance shapes in that 
range. The problem can be solved by presenting the results in narrow intervals of 
the probed phase space. Similarly, theoretical calculations of the average polarization 
in a certain experiment should consider how the momentum distribution is distorted 
by its acceptance. Alternatively, the prediction should avoid as much as possible 
kinematic integrations or provide event-level information to be embedded into Monte 
Carlo simulations of the experiments. 

The general frame-transformation relations in Eq. |2] imply the existence of an 
invariant quantity, definable in terms of A^, A^ and A^^, in one of the following 



equivalent forms: 

^ (3 + A^,) + ci(l-A^) 

^"'^ C2(3 + A^) + C3(l-A^)- ^^ 

An account of the fundamental meaning of the frame-invariance of these quantities 
can be found in Ref. [10]. We will consider here, specifically, the form 

A = ^^{-3,0,1} = ^_^ • (6) 

In the special case when the observed distribution is the superposition of n "elemen- 
tary" distributions of the kind 1 + A^ cos^i?, with event weights /*^*\ with respect to 
n different polarization axes, A represents a weighted average of the n polarizations, 
made irrespectively of the orientations of the corresponding axes: 






A-ET^^?/ErfTW P) 



The determination of an invariant quantity is immune to "extrinsic" kinematic de- 
pendencies induced by the observation perspective and is, therefore, less acceptance- 
dependent than the anisotropy parameters A^, A<^ and A^^. For instance, in the case 
illustrated in Fig. [21 as well as in the simpler case of Fig. [T]^a-b-c), any arbitrary 
choice of the experimental observation frame will always yield the value A = 1, in- 
dependently of kinematics. This particular case, where all contributing processes 
are transversely polarized, is formally equivalent to the Lam- Tung relation [TT], as 
discussed in Ref. [10]. Analogously, the example represented in Fig. [Tl^d-e-f), or any 
other case where all polarizations are longitudinal, yields A = —1. 

The existence of frame-invariant parameters also provides a useful instrument for 
experimental analyses. Checking, for example, that the same value of an invariant 
quantity (Eq. [5]) is obtained (within systematic uncertainties) in two distinct polar- 
ization frames is a non-trivial verification of the absence of unaccounted systematic 
effects. In fact, detector geometry and/or data selection constraints strongly polarize 
the reconstructed dilepton events. Background processes also affect the measured 
polarization, if not well subtracted. The spurious anisotropics induced by detector 
effects and background do not obey the frame transformation rules characteristic of 
a physical J = 1 state. If not well corrected and subtracted, these effects will dis- 
tort the shape of the measured decay distribution differently in different polarization 
frames. In particular, they will violate the frame-independent relations between the 
angular parameters. Any two physical polarization axes (defined in the rest frame of 
the meson and belonging to the production plane) may be chosen to perform these 
tests. The HX and CS frames are ideal choices at high p^, where they tend to be 
orthogonal to each other (in Eq. HJ sin^ 5 — )■ 1 for pt 3> m). At low pt, where the dif- 
ference between the two frames vanishes, any of the two and its exact orthogonal may 
be used to maximize the significance of the test. Given that A is "homogeneous" to 
the anisotropy parameters, the difference A*^^) — A^^^ between the results obtained in 
two frames provides a direct evaluation of the level of systematic errors not accounted 
in the analysis. 
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We conclude with a summary of our messages. Choosing a given polarization 
axis in an experimental analysis or theoretical calculation has a radical effect on 
the magnitudes and signs of the polar and azimuthal anisotropics observed in the 
dilepton decay distribution: all terms of the distribution must be determined to 
provide unambiguously interpretable physical information. Rotational invariance im- 
poses frame-invariant constraints on the polar and azimuthal anisotropy parameters 
and, for any mixture of production mechanisms in a given kinematic condition, there 
exists an invariant relation depending on one frame-independent parameter. Mea- 
surements of the anisotropy parameters are necessarily affected by the fact that the 
transformation from one frame to another is an explicit function of the production 
kinematics. These effects may result in a misleading interpretation of the measure- 
ments. Reporting polarization results in terms of invariant quantities facilitates the 
comparison between different measurements, and with theory, reducing the kinematic 
dependencies to their intrinsic (and physically most relevant) component. The invari- 
ant relation can also be used in the data analyses to perform self-consistency checks 
which can expose unaccounted detector effects or eventual biases. 
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